Electromechanical properties of suspended Graphene 

Nanoribbons 



Oded Hod 

School of Chemistry, The Sackler Faculty of Exact Sciences, 
Tel Aviv University, Tel Aviv 69978, Israel 

Gustavo E. Scuseria 
Department of Chemistry, Rice University, Houston, Texas 77005-1892 

(Dated: May 5, 2009) 

Abstract 

Graphene nanoribbons present diverse electronic properties ranging from semi-conductingii2j^ 
to half-metallic,— depending on their geometry, dimensions and chemical composition.— Here we 
present a route to control these properties via externally applied mechanical deformations. Us- 
ing state-of-the-art density functional theory calculations combined with classical elasticity theory 
considerations, we find a remarkable Young's modulus value of ~ 7 TPa for ultra-narrow graphene 
strips and a pronounced electromechanical response towards bending and torsional deformations. 
Given the current advances in the synthesis of nanoscale graphene derivatives, our predictions can 
be experimentally verified opening the way to the design and fabrication of miniature electrome- 
chanical sensors and devices based on ultra-narrow graphene nanoribbons. 
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The promise of graphene nanoribbons (GNRs) as ultimate building blocks for future na- 
noelectromechanical systems (NEMS) has been recently demonstrated experimentally for the 
first time.- Since their initial successful fabrication,- the dimensions of GNRs have rapidly 
reduced from the microscale down to a record breaking value of a few nanometers in width, 
fabricated by either top-down*'^ or bottom-up^S approaches. While the electromechanical 
characteristics of their cylindrical counterparts, carbon nanotubes (CNTs), have attracted 
great interest in recent years,— similar effects in GNRs have been hardly 
addressed, thus far. Here, we study the electromechanical response of GNRs by consid- 
ering a large set of hydrogen-terminated GNRs with varying lengths and widths. The effect 
of uniaxial strain on the electronic properties of infinite GNRs has been recently studied in 
detaiL^^ To simulate the bending and torsion deformations of a finite suspended GNR- we 
clamp the atoms in the region close to the zigzag edges (highlighted by orange rectangles 
in the upper panel of Fig. [1]) while applying the deformation to the atoms residing within a 
narrow strip along its center-line (highlighted by a yellow rectangle in the figure). For bend- 
ing deformations the central atomic strip is depressed down with respect to the fixed edges 
parallel to its original location in the plane of the unperturbed nanoribbon (lower left panel 
of Fig. [1]). Torsional deformations are simulated by rotating the central atomic strip around 
the main axis of the GNR and fixing its location with respect to the plane of the clamped 
edge atoms (lower right panel of Fig. [1]). The positions of the remaining atoms are relaxed 
using spin-polarized density functional theory calculations within the screened-exchange hy- 
brid approximation of Heyd, Scuseria, and Ernzerhof (HSE06), ti^^i^^ as implemented in the 
development version of the Gaussian suit of programs^--^ and the double-C polarized 6-31G** 
Gaussian basis set.— This theoretical approach has been recently shown to describe the 
physical properties of graphene based materials with exceptional success. -i^i^i^Si^i^^i^ 

To understand the mechanical response of GNRs under an externally applied stress, 
macroscopic elasticity theory concepts are here adopted. We quantify the bending deforma- 
tion by defining the depression depth, d, as the distance of the central atomic strip from the 

^The HOMO-LUMO gap energy is of central importance in the present study. Therefore, we use the HSE06 
functional, which accurately predicts bandgap values for GNRs<^i^ Values calculated using the local den- 
sity or the gradient corrected approximation usually underestimate GNRs bandgaps when compared to 
experimental measurements. 
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FIG. 1: Illustration of a graphene nanoribbon based electromechanical device. Upper panel: a 
representative structure of the graphene nanoribbons set studied. The notation MxN is used to 
represent a GNR with M hydrogen atoms passivating each zigzag edge and N hydrogen atoms 
passivating the armchair edge. The deformations are simulated by clamping the atomic strips at 
the zigzag edges (highlighted by orange rectangles) and depressing or rotating the atoms at the 
center of the ribbon (highlighted by a yellow rectangle) with respect to the fixed edges. Lower 
panels: An artist view of a suspended GNR under bending (lower left panel) and torsional (lower 
right panel) deformation induced by an external tip. Side views of the deformed structures are 
presented for clarity as insets. 



original plane of the unconstrained nanoribbon (see lower right panel of Fig. [2]). Similarly, 
the torsional angle is defined as the angle between the plane of the central atomic strip and 
the plane defined by the clamped edge atoms (see lower right panel of Fig. H]). In Fig. [2l 
the change in total energy as a function of the depression depth within the linear response 
regime is presented for the full set of GNRs studied. The marks represent calculated values 
while the lines are parabolic fits, from which the bending force constants are extracted. 
The following simple relation describes the dependence of the bending force constant of a 
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FIG. 2: Total energy changes of the 03xA^ (upper left panel), 04xA^ (upper right panel), and OSxA'^ 
(lower left panel) suspended graphene nanoribbons due to an externally applied bending stress. 
Marks represent calculated results. Lines are parabolic fits indicating that all systems are within 
the linear response regime. Lower right panel: the depression depth is defined as the distance 
between the lowered central atomic strip and the original plane of the unconstrained GNR. 

macroscopic rectangular rod on its dimensions:— 

K, = 16Yw(^^ , (1) 

where w is the width of the rod, t the thickness of the rod, / the length of the rod, and Y 
its Young's modulus. We now fit the bending force constants calculated via the parabolic 
curves in Fig. [2] to this simple relation. The length of the sample is taken to be the mini- 
mum distance between fixed edge atoms on the opposing clamped edges of the GNR. The 
width is taken as the minimum distance between hydrogen atoms passivating opposite arm- 
chair edges. For the thickness, we assume a typical value of t = 0.75 A for the graphene 
sheet^. In Fig [3] the dependence of the bending force constant on the inverse cube of the 
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length (left panel) and on the width (right panel) of the ribbon are presented. We find 
excellent correlation between our calculated values and predictions from elasticity theory. 
The Young's modulus of GNRs can be extracted from the slopes of the linear curves. An 
impressive large value of ~ 7 TPa is obtained, exceeding the measured value for micrometer 
scale suspended graphene sheets^ and the highest value calculated for CNTs with similar 
thickness parameters.— 




L"' (A"') W (A) 

FIG. 3: The dependence of the bending force constant on the inverse cube of the length at constant 
width (left panel) and of the width at constant length (right panel) of the GNRs studied. A Young's 
modulus value of ~ 7 TPa is calculated from the slopes of the linear curves by the use of Eq. [H 

A similar analysis can be performed for torsional deformations. In Fig. HJ the change in 
total energy as a function of the torsional angle within the linear response regime is pre- 
sented for the full set of GNRs studied. As before, the marks represent calculated values 
and the lines are parabolic fits. While, similar to the case of bending deformations, excellent 
parabolic fits are obtained, classical theory of elasticity fails to predict the behavior of the 
torsional force constant with the dimensions of the system and hence it is impossible to 
extract the shear modulus of GNRs from these calculations. Several factors may limit the 
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FIG. 4: Total energy changes of the 03xiV (upper left panel), 04xA^ (upper right panel), and 05xN 
(lower left panel) suspended graphene nanoribbons due to an externally applied torsional stress. 
Marks represent calculated results. Lines are parabolic fits indicating that all systems are within 
the linear response regime. Lower right panel: the torsional angle is defined as the angle between 
the rotated central atomic strip and the original plane of the unconstrained GNR. 

validity of the classical theory for the description of torsional deformations in the considered 
molecular systems. The poorly defined thickness leading to extremely large surface to vol- 
ume ratios combined with the overall nanoscopic dimensions of the systems suggest that the 
equations should be quantized in order to predict the correct behavior. Furthermore, the 
effect of elongation during the twisting process may introduce considerable deviations from 
the pure torsional equations. Therefore, it is impressive that for the case of bending defor- 
mations the effect of these factors is small and very good agreement between the behavior 
of a macroscopic rod and that of a nanoscale atomic sheet is obtained. 

Having explored the mechanical properties of suspended GNRs, we now turn to evaluate 
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their electromechanical responses. For this purpose, we consider large bending and torsional 
deformations, well beyond the linear response regime. Remarkably, we find that most of the 
studied systems can sustain extreme mechanical deformations while retaining their elastic 
nature. The HOMO-LUMO gap, which is the difference between the highest-occupied and 
lowest-unoccupied molecular orbital energies, is used to evaluate the influence of the me- 
chanical deformations on the electronic properties of the system. In Fig. [5] the energy gap 
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FIG. 5: HOMO-LUMO gap changes of the 03xA^ (upper left panel), 04x7V (upper right panel), 
and 05xA^ (lower left panel) suspended graphene nanoribbons due to an externally applied bending 
stress. 

as a function of the depression depth is presented for the full set of GNRs studied. High 
sensitivity of the energy gap upon the depression depth is found. For most of the systems 
studied, the gap energy changes by ~ 0.5 eV upon a depression of 1 nm. As expected, 
the smaller the system dimensions, the higher the sensitivity upon similar stress conditions. 
Interestingly, the 03xA^ and the narrower 04xA^ group members show an initial decrease 
in the energy gap upon bending while the 05xA^ group presents an initial increase. This 
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resembles the case of CNTs, where strain induced bandgap dependence was attributed to 
Brillouin zone deformations upon the development of mechanical stresses in the system.— 
If we regard the finite GNRs as unrolled segments of zigzag CNTs, and follow the theory of 
Yang and Han,^^ we expect to see minor bandgap changes due to bending. Nevertheless, be- 
cause of the doubly-clamped geometry, the bending is not pure and includes a considerable 
stretching component. In the case of (n, 0) zigzag CNTs, the sign of the bandgap change 
upon uniaxial stretching sgn{AEgap) = sgn{2p+l), where n = 3q+p, q being an integer 
number, and p = 0, ±1. Therefore, the bandgap of metallic CNTs {p = 0) is expected to in- 
crease upon uniaxial stretching. This is the case for the 05xA^ group, which can be regarded 
as an unrolled (6, 0) CNT and in the limit of n — oo becomes nearly metallic.-*^ When 
regarding the 04xA^ systems as unrolled (5, 0) CNTs, which are characterized by p = — 1, it 
is expected that the energy gap will decrease upon deformations which involve bond stretch- 
ing, as is the case for the 04x14 and 04x20 nanoribbons. Interestingly, deviations from this 
role occur for the 04x28 and the 03xiV group members, which correspond to an unrolled 
(4, 0) CNT with p = 1. Such deviations are expected due to the fact that the systems under 
consideration are of extremely small dimensions, where the relevance of the band-structure 
theory, upon which the Yang-Han theory is based, is limited. Furthermore, the Born- Von 
Karman boundary conditions in the circumferential direction of CNTs, taken into account 
in the Yang-Han theory are replaced by particle-in-a-box like boundary conditions in GNRs, 
which may change the overall behavior of the bandgap-strain relationship. Another interest- 
ing prediction made by the Yang-Han theory of CNTs electromechanical response is periodic 
oscillations of the bandgap with the applied strain, as the shifted Fermi points cross different 
allowed sub-bands. Evidence of such oscillations can be seen in the upper- and lower-left 
panels of Fig. [51 Due to the relatively small dimensions of the GNRs studied, only a par- 
tial period is obtained within the deformation depth range studied. This resembles recent 
measurements made on CNTs, where the frequency of the bandgap oscillations was found 
to decrease with increasing diameter.— 

A very similar picture arises for the case of torsional deformations. Fig. [6] presents the 
dependence of the energy gap on the torsional angle, up to a value of = 90°. As for the 
case of bending, energy gap changes of up to 0.5 eV are obtained upon torsion of the smaller 
systems studied. The sensitivity reduces as the dimensions of the systems are increased, and 
the general trend of increase (decrease) in the energy gap of the 03xiV, 04xA^ (05xA^), upon 
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of stresses in the system, remains. Remnants of the energy gap oscillations 
the 05x14 nanoribbon in the lower left panel of the figure. 





FIG. 6: HOMO-LUMO gap changes of the 03xA^ (upper left panel), 04xA^ (upper right panel), and 
05xA^ (lower left panel) suspended graphene nanoribbons due to an externally applied torsional 
stress. 

In summary, we have studied the electromechanical properties of suspended graphene 
nanoribbons under bending and torsional deformations. High sensitivity of the electronic 
properties to applied stresses was found, suggesting their potential use as building blocks in 
nanoelectromechanical devices. Classical elasticity theory adequately describes the mechan- 
ical behavior of the systems under bending deformations. The calculated Young's modulus 
of 7 TPa marks ultra-narrow GNRs as one of the strongest existing materials. While the 
systems we study are probably too small to be manipulated by an external tip as illustrated 
in Fig. [H the results of our calculations present important trends that are expected to 
hold for larger systems as well. An alternative platform to induce bending deformations on 
molecular graphene derivatives^ such as those studied herein, while measuring their trans- 
port properties would be the use of mechanically controllable break junctions.— In such 
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a setup a graphene ribbon bridging the gap of a pre-designed molecular scale junction may 
be subject to delicate bending deformations via the careful manipulation of the underlying 
surface. With this respect, an interesting experimental challenge would be to induce, in a 
similar manner, torsional deformations of graphene nanoribbons. 
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